Enhanced low-Reynolds-number propulsion in heterogeneous viscous environments 
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It has been known for some time that some microorganisms can swim faster in high-viscosity 
gel-forming polymer solutions. These gel-like media come to mimic highly viscous heterogeneous 
environment that these microorganisms encounter in-vivo. The qualitative explanation of this phe- 
nomena first offered by Berg and Turner [Nature (London) 278, 349 (1979)], suggests that propulsion 
enhancement is a result of flagellum pushing on quasi-rigid loose polymer network formed in some 
polymer solutions. Inspired by these observations, inertia-less propulsion in a heterogeneous vis- 
cous medium composed of sparse array of stationary obstacles embedded into an incompressible 
Newtonian liquid is considered. It is demonstrated that for prescribed propulsion gaits, includ- 
ing propagating surface distortions and rotating helical filament, the propulsion speed is enhanced 
when compared to swimming in purely viscous solvent. It is also shown that the locomotion in 
heterogenous viscous media is characterized by improved hydrodynamic efficiency. The results of 
the rigorous numerical simulation of the rotating helical filament propelled through a random sparse 
array of stationary obstructions are in close agreement with predictions of the proposed resistive 
force theory based on effective media approximation. 

PACS numbers: 47.57.-s, 47.63. Gd, 47.63. mf, 87.17.Jj 



I. INTRODUCTION 

In the past years there has been an increasing inter- 
est in propulsion on small scales, both theoretically and 
experimentally. Interest to some natural modes of low- 
Reynolds-number locomotion has a long record in applied 
mathematics [r, and underlying mechanisms of propul- 
sion powered by flexible elastic filament 1, 0) Hi j rotating 
helical flagellum 0, H, 0|, beating cilia p|, surface distor- 
tions of non-flagellated squirmers [1, [lol| and some others 
are quite well understood. Artificial nature-inspired pro- 
pellers, powered by either beating or rotating filaments 
[ll| were recently fabricated and tested vs. the theo- 
retical predictions; performance of such devices was also 
studied numerically via particle-based algorithms [T5| . 

Theoretical work on zero- Reynolds-number locomotion 
strategies (that are not necessarily biomimetic) for arti- 
ficial micro-swimmers has attracted some attention quite 
recently and several modes of propulsion, such as three- 
link Purcell's swimmer p^ . its "symmetrized" version 
[Til and generalized N-link swimmer [l^, three-sphere 
propeller [T6| . swimmer propelled by arbitrary non- 
retractable cyclic shape strokes [l3|, two-sphere "push- 
mepuUyou" |18|. surface treadmilling [T§| and surface 
tank-treading [20| and others, were proposed and stud- 
ied in details. A comprehensive review that provides the 
reader with state-of-the-art in low-Reynolds-number lo- 
comotion can be found in [21,] . 

However, the aforementioned works address the 
propulsion through an unbounded Newtonian viscous 
liquid. Some bacterial cells are used to swim though 
complex and highly heterogeneous viscous environments 



rather than Newtonian viscous liquids, such as marine 
water. For instance, clinically important spirochetes nav- 
igate efficiently through dense extracellular matrix in 
host tissues and cross the blood-brain barrier [13, [II]- 
Moreover, prospective design of artificial microrobots ca- 
pable of propulsion through soft tissues, digestion tract, 
spinal canal, etc. relies on ability to efficiently navigate 
through complex heterogeneous environment. Some ex- 
perimental observations of biological propulsion in gel- 
like polymer solutions cannot be explained in the frame- 
work of standard theories. The traditional theories for 
bacterial motion through purely viscous response predict 
that the swimming speed (at constant torque produced 
by flagellar motors) monotonically decreases with vis- 
cosity. Although the naive intuition suggests that, for 
instance, a rotating helical filament would move like a 
corkscrew when propelled through a very viscous liquid, 
this is actually never the case, no matter how large is 
the viscosity, as the zero-Reynolds-number swimming is 
purely geometric [24*]. However, the swimming speed 
of Pseudomonas aeruginosa in polyvinylpyroUidone so- 
lutions increases with viscosity up to a certain point 
and thereafter decreases 125j]. Other flagellation types 
of externally flagellated bacteria {e.g. Bacillus mega- 
terium, Escherichia coli, Serratia marcescens, Sarcina 
ureae, Spirillum serpens and Thiospirillum jenense) also 
exhibit an increase in swimming speed in viscous solu- 
tions 25|. The disagreement with traditional theories 
is more drastic for spirochetes, lacking external flagella, 
in which a helical or plane wave of the cell body moves 
backward yielding rolling of the cell body and forward 
propulsion. A remarkable feature of spirochete motility 
is that cells swim faster in a high viscosity gel-like me- 
dia than they do in low- viscosity aqueous media [2^, . 
The swimming speed of Leptospira interrogans monoton- 
ically increases with viscosity in medium supplemented 
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with methylcellulose until the viscosity exceeds 300 cP 
[26| . In low- viscosity, aqueous mediuni, T. denticola is 
observed to rotate without translating [28l |. 

Berg and Turner [2^ investigated the effect of viscosity 
on propulsion of various microorganisms and suggested 
that the propulsion enhancement in in-vitro experiments 
can be caused by loose and quasi-rigid networks formed 
by entangled linear polymer molecules such as methyl- 
cellulose. They argued that the enhanced propulsion is a 
result of flagellum pushing on the polymer network and, 
therefore, the propulsion resembles the motion of a screw 
boring through wood. Merely increasing the viscosity 
of the medium using a non-gel-forming sucrose polymer, 
such as FicoU, did not enhanced motility, which indicated 
that propulsion enhancement may be dependent on vis- 
coelastic rheology of the medium. However, recent find- 
ings of [lOl, who examined the effect of vicsoelasticity 
on propulsion for various rheological models of the liq- 
uid, showed that addition of elastic response to material 
constitutive equation does not yield the enhancement of 
propulsion. The propulsion through viscoelatsic liquid is 
always hindered (for the prescribed swimming gait), both 
velocity- and efficiency-wise, when compared to locomo- 
tion through viscous (Newtonian) solvent of the same 
viscosity. 

The first attempt to explain the propulsion enhance- 
ment theoretically, following f2§\, was made in [sij for 
externally flagellated bacteria and later extended in [32| 
to spirochetes. They suggested that the local viscous 
resistance to the motion of the flagellum normal to its 
surface is considerably higher than that corresponding 
to the tangential motion (note that in purely viscous liq- 
uid the ratio of viscous resistances of a slender body is 
^2). Indeed, the normal motion of the helical threads 
yields pushing on polymer aggregates surrounding flag- 
ellum, while tangential motions yield no perturbation to 
the microstructure. Therefore, in gel-like media with a 
sufficient traction, the flagellum is propelled much like 
a corkscrew, without slippage. On the other hand, in- 
crease in the drag on the large bacterium head results in 
decay of the propulsion speed starting at some viscosity 
in agreement with previous experiments. In accord with 
these arguments is also the monotonic increase of the 
swimming speed of the organisms lacking a large pas- 
sive head, such as spirochetes. When surrounded by a 
dense gel- like environment, the helical or flat waves of 
their body - depending on the spirochete species - get 
sufficient traction to propel the cell body through the 
medium like a corkscrew, i.e. without slippage. If trac- 
tion is not sufficient, then the cell body slips against the 
external medium and the cell translates more slowly. Al- 
though the model proposed in [3ll,[33| shows limited qual- 
itative agreement with experimental results, it involves 
two unknown phenomenological parameters - "apparent 
viscosities" of unclear physical origin. 

The aim of the present study is to develop a more rig- 
orous theoretical framework of propulsion through het- 
erogeneous viscous media. In some sense, the working 



hypothesis behind this theory follows the original Berg 
and Turner's proposal ^29J, i.e. the flagellum is pushing 
on the stationary matrix of obstacles, whereas the in- 
teraction is indirect, i.e. mediated by a viscous solvent. 
The theory is based on solutions of averaged equations 
of viscous flow through random sparse array of obstruc- 
tions such as fibers or spheres, mimicking heterogeneous 
gel-like polymer solutions where dense cores of the micro- 
gel particles (or chain-like aggregates in methylcellulose 
solutions [13) are surrounded by a viscous solvent con- 
taining dissolved polymer chains [13]. We are interested 
to address the effect of the sparse network of stationary 
obstacles {e.g. fibers or spherical cores) embedded into a 
viscous incompressible Newtonian solvent on propulsion 
of a force-(torque-)free swimmer. The properly averaged 
hydrodynamic fields {i.e. pressure and velocity) in such 
media are governed by the effective Brinkman media ap- 
proximation [ssj . 

— Vp -I- /iAti — fin ■ u , (1) 

where u is the average velocity field satisfying the conti- 
nuity equation for the incompressible liquid, V-it ~ 0, pis 
the average pressure field and k is the tensorial damping 
coefficient. The essence of this mean field model is that, 
on average, the fluid in proximity to a stationary obsta- 
cle experiences a damping body force proportional to the 
local velocity accounting for the influence of the neigh- 
boring objects on the flow. By averaged quantities we 
mean ensemble averages over all possible arrangements 
of the surrounding obstructions about a "test" obstacle. 
Far from the "test" obstacle the velocity gradients are 
weak and the equation ((!]) reduces to the differential form 
of Darcy's equation describing the flow through porous 
media and thus n is taken to be the Darcy resistance (or 
inverse of permeability). Spatially isotropic random ma- 
trix of obstacles can be described by a scalar resistance, 
Kij = Sij with 6ij being the identity tensor. Thus, 
deflnes a new length scale related to the stationary 
obstacles (usually referred as "screening" or shielding" 
length), so that for r < a^^ the velocity fleld is Stokesian, 
while for r > the velocity satisfies Darcy's equation. 
In particular, the fundamental solution of the singularly 
forced Brinkman equation {i.e. flow driven by a point 
force fi acting at Xq) are identical to these corresponding 
to (Fourier-transformed) equations of oscillatory Stokes 
flow where the frequency parameter = —'njjjv replaces 
the hydrodynamic resistance [1^ 

u,{x) = g,, {x, Xq) /, , = A ^ + B ^ , (2) 
87r/i ^ ^0 

with 

A = 2e-''(l + p-i+p-')-2p-2, 
B = -2e-'' (1 3p-i + 3p-2) + 6p-2^ 

where x = x — Xq, r = \x\ and p = ar. The funda- 
mental solution for the pressure field corresponding to 
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the Stokes flow [36[ remains unchanged. As expected, 
for a Stokes flow solution is recovered a,s A, B ^ 1. 
Thus, the flow through random matrix of obstructions 
due to a point force decays hke (due to exponen- 
tial screening of hydrodynamic interaction via e~""' in 
the above "shielded Stokeslet" solution) vs. decay 
of the classical point force (Stokeslet) solution in viscous 
incompressible liquid. 

The dumping coefflcient k in ([T]) can be deduced the- 
oretically for some representative obstacle arrangements 
by relating it to the total drag on the individual ele- 
ments comprising the network [33, . The calculation 
is based on the fact that the mean resistance force per 
unit volume is equal to (minus) the number density of 
obstructions times the mean drag force f i exerted on a 
single obstruction in the uniform flow U though a sparse 
matrix, Vp — nFi — fxa^U. For the simplest model 
system of sparse random matrix composed of stationary 
spherical obstructions of radii a, so that -Fi = dnfiaU, 
we arrive at = 90/2a^, where is the volume fraction 
or concentration of the obstacles. Self-consistent calcu- 
lations of the components of k corresponding to flbrous 
media for various spatial arrangement of long fibers, can 
be found in 
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37|. 



It should be noted that effective media approxima- 
tion ^ was originally derived in [35| from heuristic ar- 
guments. However, it was shown later to provide the 
leading correction to the Stokes drag on a test sphere 
for random distribution of stationary spherical obstacles 
in a rigorous theory where the localized resistance due 
to neighboring obstacles was taken into account in the 
multi-particle expansion scheme for small concentration 
of obstacles The excellent agreement between the 

second-order (in concentration) multi-particle theory and 
the theory based on ([T]) suggested that the latter may 
serve as an accurate approximation even in the case of 
moderate obstacle concentration. The numerical results 
based on Stokesian Dynamics suggested that ^ accu- 
rately describes the flow through the random arrays of 
spheres with up to 30% (by volume) fraction of obstruc- 
tions [soj . 

In this work we consider several locomotion modes ex- 
tensively studied in the past and relevant for propulsion 
of flagellated as well as non-flagellated organisms. First, 
we consider the G. I. Taylor's undulating deformable 
sheet propagating small-amplitude traveling waves along 
its surface Q and propelled through heterogeneous vis- 
cous medium whereas hydrodynamics is governed by 
((!]). The study of propulsion of non- flagellated swim- 
mers powered by small-amplitude surface distortions {i.e. 
squirmers) is greatly simplified via the use of Lorentz re- 
ciprocal theorem. Finally, we propose a local resistive 
force theory of propulsion through the effective media 
and test the predictions of the theory using rigorous nu- 
merical simulation of rotating helical filament through a 
sparse array of stationary spherical obstructions embed- 
ded within viscous Newtonian solvent. 



u 



c o 



X 



FIG. 1: (Color online) Schematic of the undulating sheet (red, 
solid) embedded in stationary random matrix of obstructions 
(blue dots). The sheet is propagating a traveling transverse 
wave in the positive x direction with the wave speed c = u/k 
and propelled to the left with velocity U . 



II. PROPULSION OF AN UNDULATING 
PLANAR SHEET 

A. Transverse distortions 

The 2D analog of the flagellum-powered swimming 
is the inflnite sheet along which the traveling waves 
are propagating. G. I. Taylor was the first to con- 
sider this model problem in the limit of long-wavelength 
(or, alternatively, small-amplitude) waves [3|. The same 
model was considered in 30] for modeling propulsion 
in viscoelastic liquid. The extension of the asymptotic 
small-amplitude analysis towards the case when the flow 
around the sheet is governed by (H)) is rather straightfor- 
ward. We first consider purely normal mode of deforma- 
tion, so that in the reference frame fixed with sheet the 
position of the material points follow 



y,„ = &sin [kx — Lot) , 



(3) 



where h is the amplitude of modulation, k - the wavenum- 
ber and w is the frequency of the wave. The traveling 
wave is propagating in the positive x-direction with the 
wave speed c = uj/k (the schematic of the problem is 
shown in Fig.[T]). Since the problem is 2D we recast the 
equation of motions H]) and the no-slip boundary condi- 
tions in terms of the streamfunction, if). The streamfunc- 
tion is defined by u = dyil) , v = —dxipj where u = {u, v} 
so that the incompressibility condition V • u = is triv- 
ially satisfied. Assuming that the arrangement of the ob- 
structions is random and isotropic (random arrangement 
of fibers, aligned in all directions or spherical obstacles 
in a random arrangement), the damping coefficient Kij in 
(fT|) is diagonal and equal to, a'^Sij, where a has the di- 
mensions of inverse length. Then, applying curl on both 
sides of (HI yields 



A (A - V-) = , 



(4) 



where A = V^. The no-slip boundary conditions for the 
velocity components at the deformed sheet read 







V = —buj cos (kx — Lut) 



(5) 
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Obviously, purely transverse surface distortions do not 
preserve surface area, as Vg • it ^ 0, where Vs — 
(I — nn) ■ V is the surface gradient operator, and cannot 
describe deformation of flexible but incompressible sheet. 
In the latter case, the velocity distribution is given by 
a combination of transverse and tangential distortions. 
In the leading approximation the 2D incompressibility 
condition [dus/ds + kw„ = 0, s is the arc length mea- 
sured in the direction of tangent unit vector s) yields 
dx"^ + l/'v = 0, that leads to [3[ 

u ~ —b'^koj cos {2kx — 2Lut), v — ~bu! cos {kx — cut) , (6) 

We will discuss the propulsion by tangential distortions 
separately in the next section. 

We expand the velocity at the deformed sheet in terms 
of its value at the undistorted surface (at time t = Q, 
the choice of time is arbitrary) and choose the follow- 
ing characteristic scales: k~^ for length, huj for velocity 
and {bijjk)~^ for time. In terms of the scaled variables 
we obtain the following problem for the nondimensional 
streamfunction 

A(A-a^)V = 0, y>0, (7) 
VV' + £ sinx V9j,^/; = cosx Ca; , y = 0, (8) 

where a* = a/A: is the scaled resistance and e — bk 
is the scaled amplitude of the wave. We now expand 
the solution in series of powers of e, -0 = Vo + e^i + . . .. 
Submitting this expansion into (l7][8]) and matching terms 
of the same powers e we can construct the asymptotic 
solution. The net displacement will be determined by 
the x- velocity component of the flow far from the sheet, 

Uoo = dyil> , 2/ — > oo . (9) 

In the leading order 0[e'^) produces 

A (A - al) V-o = 0, y > 0, (10) 
Vi/'o = cosa; Cj;, y = 0. (11) 

The zeroth order solution is readily found as 

= (Aoe-\/^^ + Boc-y"^ sin a; , (12) 

where the constants and Bq are determined from (jlip , 

1 



Bo = 



(13) 



The propulsion velocity, as in the case of purely viscous 
liquid appears at the next order (because of symmetry 
£ —e) as dytpo ^ as y ^ oo. 

At 0{e) we obtain the following problem 

A{A~al)^i=0, y>0, (14) 
VV'i = — sina;V9j^V'o , y = 0. (15) 



Substitution of the zeroth order solution ((T2| into ([151) 
at y = yields 

dyipi = — ain X dyyipo = 

- [{1 + al)Ao + Bo] - ^ cos 2x\ , (16) 



dxipi = - sina;9^a;'(/'o = 0, 



(17) 



as dyi/jQ = at y = 0. This suggests the form of the 
solution for ipi 



^1 =Ciy+ Mie 



Bic-^y) cos 2a; . (18) 



Substitution of ^ into (|16m7p and using (jl3p yields the 
unknown coefficients Ai, Bi and Ci, 

v/r+^(2 + 74T^) ^_ yY^ 

2ai ' ~ -^r- 

The coefficient of the linear term, Ci, gives the velocity 
far from the undulating sheet. 



Woo = dy^l\. 



--[{l + al)Ao + Bo\ = ^ 2 ■ 



Thus, the sheet is propelled to the negative x direction 
(propulsion velocity, U, is equal to minus Uoo) as in the 
case of locomotion through viscous liquid without ob- 
structions. Returning to the dimensional quantities, we 
obtain 



{bkfc 



(19) 



where c = u/k stands for the speed of traveling wave. 
The classical Taylor's results is recovered at the limit 
a ^ 0. The propulsion velocity scaled with the wave 
speed is plotted in Fig. [21 The crossover from Stokesian 
propulsion (uoo = (6fc)^c/2) to that dominated by hydro- 
dynamic resistance due to obstructions occurs at a* ~ 1. 
In other words, propulsion powered by short wavelength 
surface distortions [a/k < 1) is not affected by the ob- 
structions, while locomotion by long wavelength distor- 
tions (comparative to a~^) is enhanced, as the swimmer 
is experiencing the mean resistance due to stationary ob- 
stacles. 



B. Tangential distortions 

Let us consider now the effect of stationary obstacles 
on propulsion driven by purely tangential surface distor- 
tions (swimming gait of some non-flagellated microorgan- 
isms such as cyanobarcteria, [9i,llQi])- For this purpose we 
use the same model problem of the 2D sheet propagat- 
ing traveling waves of compression/extension along its 
surface in the positive x direction. 



Vn 



y ■ 



Xm = X + b sin {kx — ut) 



(20) 
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and the remaining constants are readily found 



= -Bi = - 



1 
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The scaled fluid velocity at y oo is the same 
as for purely viscous liquid with viscosity Uoo = 
dy'>l'i\{x,y^oc) — ~£/2, and being re-written in a dimen- 
sional form yields 



c{hkf 



(25) 



FIG. 2: Propulsion velocity Uac/{bk)^c vs. the scaled mean re- 
sistance due to stationary obstructions, a/k (log-linear plot). 
The crossover from Stokesian regime to the enhanced regime 
dominated by the network hydrodynamic resistance occurs 
when screening length becomes comparable to the wave- 
length k^^ . 



so the displacement is in the direction of wave propa- 
gation (positive X direction). Therefore, the propulsion 
of the planer sheet powered by purely tangential small- 
amplitude surface distortions is not affected by presence 
of obstacles in the first approximation. 



Similarly to the analysis of the transverse waves (Eqs. 
[7][51) , for purely tangential traveling waves we obtain the 
following problem in terms of the dimensionless stream 
function ijj at the undeformed surface: 

A(A-a2)^ = 0, 2/>0, (21) 
Vip + s sinx \7dxip = — cos X By, y — 0. (22) 

Substituting the asymptotic expansion ip — 'ipo + sipi + 
. . . into (|21j|22p results at the leading order 0(e°) in the 
following boundary conditions 



Vi/jo — — cos X By , y — . 



while at any order -!/;„ should obey (|2T|) . 
The zeroth order solution is 



(23) 



(24) 



where the constants are determined from by (|23p and 
read 



Ao = -Bo = 



1 



The net propulsion velocity is again restricted to the next 
order 0(e) as dytpo ^ as y ^ oo. 

At 0{e) the boundary conditions for -01 at y = read 

V-ipi = — smxVdxipo , y = . 
Substituting the coefficients Aq, Bq we find that at y = 

dytpi = -(v/TTof^o + So) sin^x = - Q - i cos2a;^ 
dx'ipi — (Aq + i?o) cosa; sina; = 0. 

These boundary conditions suggest the form of the solu- 
tion for tpi satisfying Eq. [2T] 



C. Rate of viscous dissipation 

Since normal distortions are expected to facilitate 
propulsion, it is instructive to determine if such augmen- 
tation is costly in terms of power invested in swimming. 
The power expended by an arbitrary shaped organism 
and dissipated by viscosity is 



V = - (<T 

'S 



n) ■ u dS 



(26) 



If decompose the surface velocity as u = u' + U where 
U is the propulsion velocity of the swimmer and u' is 
the velocity of the swimming stroke, and since a ^ cr' 
the integral in (|26p remains unchanged when written in 
terms of it' as the swimmer is force-free, jg{cr-n)dS — 0, 
(we omit the prime hereafter for simplicity). Re- writing 
the dissipation integral as P = 2iJ, Jy E : E dV, where 
E is the rate-of-strain tensor and expressing the product 
E:Ea.sJ2 + '^{diUj){djUi), where ^ = curlit denotes 
vorticity, allows expressing V as [l0| 



C^dV 



2/1 / Tii {uj djUi) dS . 



(27) 



Here V is the fluid volume surrounding the swimmer. 

For purely tangential surface distortions the second 
term on r.h.s of (|77|) can be re- written as 2/i v^Ks dS, 
where = —{ds/ds)- n is the curvature of the surface 
along the direction of the flow [10], and it is zero for 
un undistorted planar sheet. For purely transverse sur- 
face distortions the surface integral on the r.h.s. of ([77|) 
is equal, to the flrst approximation, to 2/i U dxV dS, 
where U ~ ce^ and dxV = 0{cek) and, therefore, its con- 
tribution to the dissipation rate is restricted to 0{e^iic^). 
Thus, for both swimming gaits (normal and tangential 
surface waves), the rate-of-work, to the leading approx- 
imation, is given by the volume integral in l|27p over 
squared vorticity, = (A-f/j)^, 



(28) 
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Substituting ipQ corresponding to purely normal (|12p and 
purely tangential gaits, respectively, and integrating 
over volume we find the rate of viscous dissipation to the 
leading approximation. It appears that to the leading 
approximation (with an error of 0{e)) the dissipation 
rate (per unit area of the sheet) is the same for both 
swimming gaits and equal to 



2^1 



(29) 



where e — bk. The ratio of rate-of-work (per unit area 
of the sheet) in propulsion through matrix of obstruc- 
tions and that in unbounded viscous liquid {Vs — ^ic^e^k 
d, subscript 'S' stands for "Stokesian") 'P/'Ps>l and 
it is a monotonically increasing function of q,. Again, 
the effect of the embedded obstruction matrix on rate of 
viscous dissipation is evident for distortions with wave- 
lengths (fc~^) comparable to the screening length . 

The routine definition of swimming efficiency based 
on the ratio between the work invested in dragging 
the immobile swimmer and that expanded in swimming 
(Lighthill's efficiency |JJ) is not applicable for 2D swim- 
mers as the drag force is not defined due to Stokes para- 
dox. Thjs may be considered as a mere issue of normal- 
ization [19j, [2(| and a natural measure for the propulsion 
efficiency of an undulating sheet can be given by the ratio 

, nU'^k 

Using the expression (|19p for the propulsion velocity by 
the normal distortion and (1291) we arrive at 



6 = 



(1 



1,2^3/2 



2^1 



(30) 



The propulsion efficiency S/e^ in ([50]) is depicted in 
Fig. [3] as a function of a* . It is readily seen that the 
propulsion through matrix of obstructions is more ef- 
ficient that swimming in the unbounded viscous liquid 
with ^5 = We arrived at the improved propul- 

sion efficiency for transverse distortions and retarded ef- 
ficiency with longitudinal distortions. Therefore, propul- 
sion through heterogenous viscous environment can be 
advantageous both speed-wise and efficiency-wise. Al- 
though improved swimming speed is expected as a re- 
sult of hydrodynamic interactions with obstacles for a 
prescribed swimming gait, the improved propulsion effi- 
ciency is rather surprising result. A similar occurrence of 
enhancement of dragging efficiency due to hydrodynamic 
interaction between a passive load and a micro-swimmer 
towing it was found in (4(| . 



III. PROPULSION OF NON-FLAGELLATED 
SQUIRMERS BY SURFACE DISTORTIONS 

We expect similar propulsion augmentation for non- 
flagellated swimmers (such as cyanobacteria ) mov- 
ing through the heterogeneous viscous media. In this 
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FIG. 3: Propulsion efficiency of the oscillating sheet, 5/e vs. 
the dimensionless resistance, a/k (log- log plot): transverse 
surface waves (solid line); tangential (extension/compression) 
surface waves (dashed line). 



case the analysis can be greatly simplified via the use of 
Lorentz reciprocity ^36 ], that can be shown to hold for 
the effective media equations ([T]) linear in u: for any two 
arbitrary solutions (m, cr) and (m, ct) decaying far from 
the swimmer: 



/ {a ■ n) ■ u dS — \ {<j ■ n) ■ 

Js Js 



udS. 



(31) 



Here S is the instantaneous surface of the swimmer, cr ■ n 
is the local drag force the fluid exerts on S. The aim of 
this section is to show that the propulsion of cyanobac- 
teria through viscous liquid with embedded network of 
obstructions can be enhanced when locomotion is pow- 
ered by traveling normal distortions of the outer surface. 
Following [Tol| we will consider asymptotically tractable 
case of nearly spherical squirmer of radius a propelled by 
small amplitude surface waves. Substituting u = U + u', 
where U is the propulsion velocity and u' is the sur- 
face distortion in the frame of reference fixed with the 
squirmer's center-of-mass, into ([5Tt and setting the net 
force on the swimmer to zero, we arrive at 



F U 



- [ (? • n) • m' dS* , 
Js 



(32) 



where a is the stress field corresponding to the trans- 
lation of the spherical squirmer when acted upon by an 
external force F. The equality in ([311) is exact when 
the propulsion is governed by purely tangential distor- 
tions. For an auxiliary problem of towing the immobile 
spherical squirmer through viscous liquid with embedded 
sparse network of obstructions, we use the well-known re- 
sults for the drag force 



F = 6TTfj.a {l + aa+ {aaf/9) U 



and 



cr • n\ 



2a 



(1 + aa)I -\- nn 



for the local traction on the immobile squirmer's sur- 
face {e.g. [1^). The resulting equation for the propulsion 
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FIG. 4; Propulsion velocity of a spherical squirmer of radius a 
moving through heterogeneous viscous medium as compared 
to Stokesian velocity of propulsion in unbounded viscous liq- 
uid, U/Us, vs. the dimensionless resistance, aa: purely nor- 
mal surface distortions (solid line); purely tangential surface 
distortions (dashed line) 



speed in terms of surface motions reads 



U 



1 



47ra2 [1 -I- aa -f (aa)2/9] 



(1 + aa)u' + {u' ■ n) n 



AS. (33) 



The propulsion speed of a spherical squirmer in an un- 
bounded viscous liquid is recovered in the limit aa — > in 
accord with [To^, Us ~ — (47ra^)~^ u'dS (the subscript 
"S" stands for "Stokesian" ) . Decomposing the arbitrary 
surface velocity into the normal and the tangential com- 
ponent. It' = (J — nn) ■ u' + {u' ■ n)n, allows Eq. [33] to 
be re-writtcn for cither purely normal (superscript n) or 
purely tangential (superscript s) surface distortions as 



U"^' ^T'^'^iaa) Us 



(34) 



where T^^-^ = {i+-+(-)V3, i+aa} . clearly, 1< .F" < 3 
and < JF" < 1 for all values of aa > 0; the compari- 
son of the propulsion velocity in a heterogeneous viscous 
medium and a purely viscous solvent for both kinds of 
surface motion (normal and tangential) is depicted in 
Fig. m Thus, when both swimming gaits (normal and 
tangential waves) of a squirmer propelled through hetero- 
geneous viscous environment are compared to propulsion 
through unbounded viscous liquid, we find that normal 
surface distortions enhance the locomotion, while motion 
powered by purely tangential waves is hindered. 



IV. PROPULSION OF A ROTATING HELICAL 
FILAMENT 

A. Modified resistive force theory for propulsion 
tlirougli lieterogeneous viscous media 

In this section we develop a modification of the resistive 
force theory (RFT) of the propulsion through effective 



heterogeneous viscous media powered by rotating helical 
flagellum. We assume again that small inertia-less flow 
in such media is governed by ^ and we aim to find the 
propulsion velocity and rate-of-work of the rotating rigid 
helix as a function of its geometry {i.e. pitch angle and 
the radius r, see Fig. [S]) and angular velocity lu. 

The result for the propulsion speed of the force-free 
[45| slender helix rotating in Newtonian incompressible 
viscous liquid reads {e.g. 0, EH) 



U 



sin 20 



2 (1 + sin^ 



(35) 



so that the ratio of velocities of a helix to a corkscrew is 
independent of viscosity and by ([55)1 : U/ {ruj cot 6) < 1/2. 
Thus, in the unbounded Newtonian viscous liquid the 
helix needs at least two turns to progress over a distance 
of its threads. The rate-of-viscous dissipation is given by 



1 



sm 



(36) 



where f± and /|[ the force densities corresponding to 
transverse and axial motion of a filament, respectively; 
f± w 2/|| for purely viscous liquid. 

Wc arc interested to extend these results and derive 
the propulsion velocity and power as a function of the hy- 
drodynamic resistance of the random matrix of obstacles 
a. Let us consider propulsion speed of a single rigid helix 
that rotates with some constant angular speed uj. Taking 
63 be the direction of the helical axis, the helix centerline 
is given by r(s, t) = {r cos {ks + uit), r sin {ks + cot), bs + 
Ut}, where k ~ 2tt/£ with i being a length of a single 
helical turn so kr = sin 9; b ~ cos 9 and U is yet undeter- 
mined propulsion velocity. The basic assumption of the 
local RFT, is that the local force per unit length (force 
density) exerted on the slender filament is given by 



/ = f±{u- {U -3)3)+ • S)S , 



(37) 



where s = dr/ds is the local tangent, u = dr/dt is 
the local velocity and f± and /|| the force densities 
corresponding to transverse and axial translation of a 
straight cylinder. For purely viscous liquid, f± = 2/|| = 
AiTfiE, with a small parameter E = (ln2/e)~^ where 
e=:a/6<Clisan aspect ratio {2b is the cylinder length 
and 2a is its diameter). The slenderness of the fila- 
ment is controlled by the parameter Kfa <C 1, where 
Kf = \d'^r/ds^\ — SIT? 9/r is the local curvature of the 
filament centerline. 

We substitute the local velocity as it = u± -f Uej, into 
([57)) and setting the total force in the Xa-direction to zero, 
we find that 

U [/± + (/|| -/±(s-e3)')] ds = 
Jo 

-if\\-f±) I {3 ■ u^){s ■ e^)ds . (38) 
Jo 
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Substituting ^ = f±/f\\ and the expressions for the local 
tangent s and u± into (j38p one can derive the propulsion 
velocity of the force- free helical swimmer rotating around 
its central axis with angular velocity w as a function of 



U _ 1) sin 261 

~ 2[l + (C-l)sin2 e] 



(39) 



The RFT result ((35|) is readily recovered from ((39|) for 
/x = 2/|| and ^ = 2. Similarly, the rate-of-dissipation 

V = — f ■ uds, can be found 



V 



2/_lL r2w2 



l + C + cos 26* 



(40) 



that reduces to the known RFT result ([36]) for ^ = 2. 
Also, the propulsion efRciency can be defined as a ratio of 
power required to drag the helical filament along the x^- 
axis with velocity U and the power expanded in force-free 
swimming, S = F ■ U/P. The force required to tow the 
helix is easily obtained by integrating the local force / in 
dlZl) for M = C/eg: F = U {f±{l-cos^ 9) + cos^e)Le3. 
Then, using F and V yields an expression for propulsion 
efficiency as a function of the pitch angle and the ratio ^: 



- 1)2 sin^ 29 
4C 



(41) 



For unbounded viscous liquid (, = 2 and S ^ ^ sin^ 29, 
yielding the optimal efRciency of 12.5% at the pitch angle 
of 6* = 45°. 

Now, let us solve the analogous problem for the force- 
free rotating helix propelled through a heterogeneous vis- 
cous medium. In this case (1551) still holds, while the force 



FIG. 6; Propulsion of a rotating helical filament through 
heterogeneous viscous media: (a) optimal velocity, U/ur, 
vs. the scaled resistance, aa; (b) optimal pitch angle 9 — 



: arccos 



the scaled resistance aa. 



densities fj_ and /| | for a rigid cylinder of radius a trans- 
lating through medium with the effective resistance a are 
now functions of aa 37], 



- — = - [aa) - 

4:TTfi 4 



Ki{aa) 
Ko{aa) 



,/|| 1 Ki{aa) 
— — aa ——— r , 

47r^ 2 KQ{aa) 

(42) 

where a is the radius of the filament and Kp(x) are the 
modified Bessel functions of degree p. To determine the 
propulsion speed we only need to know the ratio of the 
resistance coefficients. 



CM Kojaa) ^ 
2 ' ^ - 



K\ (aa) 



(43) 



Note that in the limit of vanishing matrix resistance 
aa — > the correct viscous limit, ^ = 2, is recovered [46l |. 
Substituting ^(aa) into ([55)1 yields the expression for the 
propulsion velocity of slender helix rotating in viscous 
liquid with an embedded sparse matrix of obstacles. For 
vanishing resistance, ^ — > 2, and the propulsion velocity 
in ([55]) tends to that in ([55)) . as expected. In the opposite 
limit of large values of ^ the motion resembles that of a 
corkscrew, U = ur cot 9, boring through solid without 
slip. Optimizing swimming speed for prescribed rota- 
tion velocity, we find that swimming speed is maximized 



at the pitch angle 9 = ^ arccos 



vanishing resis- 
~ 35.26° 



tance yields the pitch angle 9 — arccos (-ij 
that maximizes the propulsion speed in purely viscous 
liquid, giving U/ujr ~ 0.354 in agreement with previous 
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FIG. 7: (Color online) Propulsion efBciency S (|41|) of a ro- 
tating helix in heterogeneous viscous media as a function of 
pitch angle 9 (rad) and hydrodynamic resistance aa. 

theories {e.g. see Q). The optimal propulsion velocity 
of the rotating helix moving through heterogeneous vis- 
cous environment is depicted in Fig. [6la as a function of 
dimcnsionlcss hydrodynamic resistance aa. The optimal 
pitch angle, however, decreases with the increase in aa 
(see Fig. |6?)). 

Lastly, the efficiency of propulsion can be found by 
substituting <^ in (j43|) into (I4ip . The resulting depen- 
dence on the scaled hydrodynamic resistance aa and the 
pitch angle 9 is depicted in Fig. [71 It can be readily 
seen that, again, the propulsion is advantageous not only 
speed-wise, but also efficiency-wise as S grows with aa. 
Interestingly, the optimal pitch angle remains 45° for all 
values of aa since S oc sin^ 29. The enhanced efficiency 
predicted by the theory is not sufficient, however, for 
maintaining the swimming speed for fixed torque applied 
to the filament. For a prescribed power expended in 
swimming, the swimming speed decays with aa, which 
indicates that in addition to the proposed enhancement 
due to embedded matrix of obstacles, active control of 
motility by microorganisms may be involved. 



B. Numerical simulations of propulsion through 
heterogeneous viscous medium 

To validate the theoretical predictions of the preced- 
ing section we implement the numerical scheme based 
on multipole expansion of the Lamb's spherical har- 
monic solution of Stokes equations (4^ . The filament is 
constructed from nearly touching rigid spheres ("shish- 
kebab" model) and the obstruction matrix is modeled as 
random sparse array of stationary spheres (see Fig. [5]). 
All N spheres (composing the filament and the matrix) 
have the same radius a. The no-slip condition at the 
surface of all spheres is enforced rigorously via the use 
of direct transformation between solid spherical harmon- 
ics centered at origins of different spheres. The method 
yields a system ofO{NC-^) linear equations for the expan- 
sion coefficients and the accuracy of calculations is con- 
trolled by the number of spherical harmonics (i.e. trunca- 
tion level), £, retained in the series. The same approach 
was used in [23, [4^ for modeling Purcell's toroidal swim- 




FIG. 8: (Color online) Illustration of a "shish-kebab" filament 
propelled through heterogeneous viscous medium modeled as 
sparse matrix of stationary spheres immersed into Newtonian 
viscous liquid. 



mer. 

The spheres composing the helical filament are par- 
titioned along the backbone of the filament, r(s) = 
{rcos(fcs), rsin(fcs), bs}, so that the distance between 
centers of neighboring spheres is set to 2.02a. Here 
kr = sin 9 and b = cos 9 with r and 9 being a radius 
and a pitch angle of the helix, respectively. The motion 
of the ith sphere composing a rotating helix can be de- 
composed into translation and rotation about its center 
as Vi = Ui + f2i X Ti with Ui = Ue^ + uie^ x Ri and 
Hi — Loe^] here Ri is a position vector to the ith sphere 
center in the fixed laboratory frame and is the radius 
vector with origin at the center of ith sphere, uj is the an- 
gular velocity of rotation and U is the propulsion velocity 
along xa-axis. The velocity of propulsion is determined 
by setting the net force exerted on rotating helix in the 
direction of propulsion to zero, Fi^ = 0, while trans- 
lation and rotational velocities of the stationary obstacles 
are both set to zero. The rotation is powered by an ex- 
ternal torque. The rate-of-work expended in propulsion 
of a rotating filament is found from 

i=l 

-^a;(-i?,;2F,;i+i?,iF,2+T,3), (44) 
1=1 

where Fi — Jg^ cr-ndS and Ti — Jg^ ri x {cr-n)dS are 
the hydrodynamic force and torque, respectively, exerted 
on ith sphere composing the filament. 

We initially test the scheme for the case of rotating 
helix moving through unbounded viscous liquid without 
obstacles. The results of the calculation are presented 
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e (rad) 

FIG. 9: (Color online) Propulsion velocity U/cor of a rotat- 
ing helical filament built of 30 spheres vs. the pitch angle 9. 
The void symbols represent propulsion in unbounded viscous 
liquid: r = 2a (□), r = 3a (Q) and r ^ Aa (A); the fuU 
symbols stand for propulsion through viscous liquid with an 
embedded random matrix of 8% (vol) stationary spherical ob- 
stacles: r = 2a (■), r = 3a (•) and r — Aa (a). The symbols 
stand for the mean value based on 20 random configuration 
of obstacles; the error bars length is doubled mean standard 
deviation. The solid line is the theoretical result for propul- 
sion in unbounded viscous liquid (|39p (with ^ = 1.52 as found 
from simulations for a "shish-kebab" rod with Np — 30, see 
Fig. I10|l and the dashed curve is the theoretical prediction 
P9|) for propulsion through sparse random matrix of obstruc- 
tions with aa — 0.6 (corresponding to a random matrix of 
spherical obstacles with <j!> — 0.08) for (,{aa) given in Eq. |43] 
and corrected for the "shish-kebab" shape of the filament. 



in Fig. [5] for a helix composed of 30 spheres with radii 
r — 2a, 3a and 4a upon varying the pitch angle 9 (empty 
symbols). The results of the calculation of the scaled 
propulsion velocity U/wr vs. 9 for small pitch angles 
(i.e. for Kftt — (a/r) sin^ 9 <C 1, where kj is a lo- 
cal curvature of the filament centerline) are in a good 
agreement with the RFT result (solid line) corrected for 
a "shish-kebab" shape of the helix. Actually, for such 
"shish-kebab" filament variation of the shape occurs on 
a scale of the filament radius and not the length and 
therefore one should not expect the result of the local 
theory based on assumption of gentle variations of the 
shape, to hold in this case. However, this only changes 
the numerical coefficients in front of the resistance co- 
efficients f± and /||, leading to a slightly different ra- 
tio f±/f\\ ~ 1.67 at e ^ (instead of ~ 2 for slender 
particles with gradual variation of the geometry). The 
coefficient can be found numerically by computing the 
two force components, transverse and longitudinal, ex- 
erted on a "shish-kebab" rod upon varying its length [i.e. 
number of spheres composing the filament). The results 
of the calculation are presented in Fig. \W\ together with 
the best fit of the form {ci + C2e) / {c^ + C4e) . The theoret- 
ical prediction of the helix velocity in Fig. [5] (solid curve) 
is based on ((39)) with ^ = f±/ f\\ = 1.52 corresponding to 
a straight "shish-kebab" filament composed of No = 30 
spheres (see Fig. [TU|. 

It can be readily seen that the agreement between the 
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FIG. 10: Numerically calculated force ratio f±/f\\ for a 
"shish-kebab" straight rod of length 26 made of A'^ spheres 
of radii a as a function of the aspect ratio e — a/b (void sym- 
bols). The solid line stands for the best fit {ci + C2e) / {c3+C4f-) . 



numerical (void symbols) and theoretical result is very 
good for pitch angles 9 < 0.3 rad, as the data correspond- 
ing to different helix radii (r = 2a, 3a and 4a) collapse 
on the theoretical curve (solid line in Fig.lHl). For higher 
pitch angle the assumption of slenderness (k/c <C 1) is 
violated and deviation from the theory is evident due 
to hydrodynamic interaction of the helix threads. In- 
crease in the radius of the helix (for the same pitch an- 
gle) yields better agreement with the RFT prediction 
for both: viscous liquid and heterogeneous viscous me- 
dia. The position of the maximum is observed for the 
pitch angle 9 ~ 0.6 (rad) in agreement with theoretical 
predictions regardless of the total length and radius of 
the helix. Increasing the number of spheres composing 
the helix results in slightly increased propulsion velocity, 
while the increase in accuracy level, £, slightly diminishes 
the speed. Note that for slender helices the calculations 
are rather accurate even for just two harmonics (£ — 2) 
retained in the series, as the hydrodynamic interaction 
between treads is negligible. 

Next we consider propulsion through a sparse (8% by 
vol) matrix of No = 30 stationary spherical obstruc- 
tions. Nonintersecting stationary spheres are arranged 
randomly in a rectangular box with its longer size ori- 
ented along the axis of the filament (see Fig. [5]). The size 
of the box is determined by the matrix density, i.e. vol- 
ume fraction (p of the obstacles in the box and the slender- 
ness of the filament. It should be realized, however, that 
the filament is propelled through a finite heterogeneous 
region surrounded by unbounded viscous liquid. How- 
ever, the hydrodynamic disturbance from the rotating 
force-free filament in the heterogeneous effective media is 
expected to decay considerably faster due to "shielded" 
interaction when compared to unbounded viscous liquid 
{u = o(l/r^), see Eq. [5]), and thus the effect of the vis- 
cous domain outside the box is expected to be negligible. 
The propulsion velocity is determined from static rather 
than dynamic calculations in the same way as for fila- 
ment rotating in purely viscous liquid. The value of the 



11 




0.10 ' ' 

0.00 0.02 0.04 0.06 0.08 

<P 

FIG. 11: Propulsion velocity of a rotating helical 'shish- 
kebab" filament (as in FigHJ with r ^ 2a, 6 = 0.3 (□) and 
r = 3a, 9 = 0.5 rad (a) vs. concentration (j> of stationary 
spherical obstructions, 0. The solid curves are the modified 
RFT result ((39)1 with i{aa) given by (|43)) . Error bars show 
the doubled mean standard deviation of 20 random configu- 
rations. 



velocity is averaged over 20 independent configurations 
and the results (full symbols) are shown in Fig. [9l It can 
be readily seen that swimming through heterogeneous 
domain (see Fig. [5]) yields faster propulsion when com- 
pared to swimming in purely viscous liquid (for the pre- 
scribed swimming gait, i.e. lo). The agreement between 
the modified RFT (dashed curve) and the results of nu- 
merical simulations (full symbols) is very good. The the- 
oretical prediction is based on ([M)) with ^ given by ([33]) 
for aa = 0.6, corresponding to obstacle concentration of 
4> = 0.08 (aa = S-^/^/^ for random array of spherical ob- 
stacles [38[), and multiplied by a factor of 0.83 correcting 
for the "shish- kebab" shape of the filament. 

The concentration dependence on the propulsion speed 
is depicted in Fig. [TT] for a helix composed of = 30 
spheres with r = 2a, = 0.3 and r — 3a, — 0.5 for con- 
centration of spherical obstacles in a cell up to 8% (vol). 
The monotonic increase in the propulsion speed U /ojr is 
evident and the agreement with the theoretical predic- 
tion ([55)1 is very good. The deviation from the theoreti- 
cal prediction at small is probably due to the fact that 
derivation in [sj leading to (|43)) breaks down when the 
screening length and the filament length become compa- 
rable; this case requires a special consideration and will 
be addressed elsewhere. 

The power, V , expended in propulsion of the "shish- 
kebab" filament is calculated via ((44|) . The correspond- 
ing hydrodynamic efficiency, (5, is determined as 5 = 
TZfuU^ I'P , where TZfu is the appropriate hydrodynamic 
resistance equal to the drag force on the (non-rotating) 
helix dragged along xa-axis with U — \. The comparison 
of the numerical results for a helix composed of 50 spheres 
and propelled through an unbounded viscous liquid, with 
the prediction of the RFT (|^T|) corrected for the "shish- 
kebab" shape [i.e. = 1.56 corresponding to Np — 50) 
is depicted in Fig. [T^] for several helix radii. We see that 
the theoretical prediction is less accurate than that for 
the propulsion speed (see Fig.lH]), and that the hydrody- 
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FIG. 12: (Color online) Efficiency of propulsion, S, of a ro- 
tating helical filament propelled through unbounded viscous 
liquid vs. the pitch angle 9. Symbols stand for the results of 
the numerical calculation for a "shish-kebab" filament com- 
posed of 50 spheres with: r = 6a (□), r = 8a (o) and r = 10a 
(A). Full squares (■) are the results for a shorter helix with 
Np — 30 and r = 6a. The continuous curve is a prediction 
of the RFT theory (I41|) with ^ — 1.56 corresponding to the 
"shish-kebab" helix composed of Np = 50 particles. 



namic interaction between treads at small r diminishes 
the propulsion efficiency when compared to the theory. 
For r = 10 a the agreement is good for pitch angles < 0.4 
rad. The position of the peak efficiency is slightly below 
the theoretical prediction of ^ « 0.79 rad and is around 
0.7 rad. Again, the increase in the radius of the helix 
(keeping helix radius fixed) results in a better agreement 
between the simulation results and the RFT prediction 
both for purely viscous liquid and heterogeneous media, 
as the slenderness parameter K/a diminishes. The effect 
of the total length of the helix on propulsion efficiency is 
quite minor, a result for helix composed of 30 spheres is 
very close to that of a longer helix with Np — 50 (see full 
squares vs. empty squares in Fig. I12p for a wide range of 
the pitch angles. 

Next we consider propulsion efficiency through sparse 
random array of stationary spherical obstructions im- 
mersed into incompressible viscous liquid. Here we use 
helices composed of 30 spheres and obstruction matrix 
composed of No = 30 spheres with volume fraction of 
obstacles of 5% (vol). The results for the swimming effi- 
ciency are shown in Fig. [13] as full symbols for r = 4a(B) 
and r = 6a (•). In accord with the prediction of the 
modified RFT (the dashed line), the helix is a more effi- 
cient propeller in the presence of obstacles: the maximum 
efficiency increases more than two folds. As for propul- 
sion through purely viscous liquid (void symbols vs. solid 
line), the theoretical prediction based on effective media 
approximation overestimates the efficiency for large pitch 
angles, where hydrodynamic interaction of the treads is 
important, while the agreement is reasonably close for 
pitch angles 9 < 0.4 rad. The RFT prediction was cal- 
culated from ([41]) using aa = 0.474 corresponding to ob- 
stacles' volume fraction oi 4> — 0.05, and ^ multiplied by 
a constant factor of 0.83 correcting for the "shish-kebab" 
shape of the filament in the same way as done previously. 
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FIG. 13: (Color online) Efficiency of propulsion, S, of a ro- 
tating helical filament composed of 30 spheres and propelled 
through viscous liquid with an embedded random matrix of 
5% (vol) stationary spherical obstacles, with vs. the pitch 
angle 9. Full symbols stand for numerical results for propul- 
sion heterogeneous viscous media with 5% (vol) of spherical 
obstacles embedded into viscous solvent: r = 4a (■), r = 6a 
(•); void symbols are the results of the numerical calculation 
for same filament, propelled through unbounded viscous liq- 
uid: r — (□), r = 6a (o). Error bars depict the doubled 
mean standard deviation based on 20 random configurations 
of the obstruction matrix. The curves stand for predictions 
of the RET theory using (|4ip for the propulsion efficiency 
in both cases: unbounded incompressible viscous liquid for 
^ = 1.52, corresponding to a helix composed of Np = 30 
spheres (solid curve) ; heterogeneous viscous media using ^ as 
in (|43ll with aa = 0.474 corresponding to obstruction volume 
fraction cf) = 0.05 (dashed curve). 



It should be stressed that no adjustable parameters are 
involved in comparison of the results of numerical simu- 
lations and the theory in Figs. [H [TTl [12] and [13] (besides 
the shape factor for the "shish-kebab" filament that is 
estimated numerically). 

Lastly, we calculate the propulsion velocity of a helix 
scaled with the rotation speed controlled by the power in- 
vested in swimming (i.e. torque applied to the filament), 
given by ^/pJJIL. The modified RFT theory (|39l40p in- 
dicates that for a p rescribed power the scaled propulsion 
speed, U / y/V/fiL, should monotonically decay with the 
increase in aa, while at aa — *■ 0, this velocity grows un- 
bounded as /|| in (I42p vanishes (while their ratio ^ 

remains finite). Note that U/ y/V/JiL (in comparison to 
U/usr and 5 in ((39|) and (|4ip. respectively) is no longer 
a sole function of the ratio ^ — f±/ f\\, but also depends 
on the force density f± via V. Obviously, as concentra- 
tion of the obstacles tends to zero, it is expected that the 
force density f± should tend to a finite Stokesian value 
corresponding to an unbounded viscous liquid. This un- 
physical behavior near aa = can be corrected by con- 
structing a proper slender body approximation for the 
Brinkman equation (|T]) and will be addressed elsewhere. 
However, the growth of U/ ^JV / jiL as aa diminishes in- 
dicates the propulsion speed can be enhanced for a fixed 
power. For this matter, we compute the scaled propul- 
sion speed of the rotating helix composed of 30 particles 
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EIG. 14: Dimensionless propulsion velocity, U / y^V/fiL, of a 
rotating "shish-kebab" helix of radius r = 4a and A'^ — 30, 
vs. the pitch angle 0. Full symbols (□) stand for a helix 
propelled through a unbounded viscous liquid of viscosity fi 
and empty symbols, (■), correspond to the helix propelled 
through heterogenous medium modeled as 5% (vol) random 
array of stationary spherical obstructions of radius a embed- 
ded into a viscous solvent of viscosity fi. The solid and the 
dashed lines stand for the RFT result of propulsion in purely 
viscous solvent and viscous liquid with an embedded matrix 
of obstacles, respectively. 

and the embedded matrix composed of No — 30 spheres 
with the volume fraction of obstacles of 5%. The re- 
sults of the computation are provided in Fig. [14] vs. the 
pitch angle, 6. The numerical results (symbols) demon- 
strate the propulsion augmentation (up to ~ 18 % of the 
Stokesian velocity) for rotating helix with fixed rate-of- 
work. The continuous lines refer to the modified RFT 
theory corrected for the shape of the filament as before, 
and the agreement with the numerical results is quite 
good for sm all pitch angles. Since the scaled propulsion 
U / yV/jlL decays with aa and in the limit of vanishing 
resistance should yield a Stokesian result, it is expected 
that there should be the optimal resistance aa maximiz- 
ing the propulsion speed of the helix. 



V. CONCLUDING REMARKS 

We have demonstrated that transverse 
strokes/motions can lead to enhancement of propulsion 
via the heterogeneous viscous environment, modeled 
as sparse matrix of stationary obstacles embedded into 
incompressible viscous liquid, for different swimming 
gaits. The tangential mot ions, however, do not offer such 
enhancement. The present theory is in accord with 
the qualitative argument made in [29j suggesting that 
"...efficiency of flagellar propulsion is probably enhanced 
by the gel-like structure" of the medium. Combining 
the local theory and rigorous numerical calculations we 
demonstrated that the propulsion can be enhanced both 
speed- and efficiency-wise. For instance, the translation 
speed of the rotating helical filament is increased not 
only for a prescribed rotation velocity but also for a 
fixed torque applied to the filament. These findings can 
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provide the physical grounds for explaining experiments 
showing enhanced propulsion of externally flagellated 
bacteria [1^, [2§| in viscous polymer gels. 

We should next address the question of applicability of 
the effective media approximation to biologically relevant 
cases. Originally, the theory was developed to described 
the resistance to the flow through such media, i.e. sta- 
tionary arrays of obstructions (fibers or spheres), while 
we consider a propulsion of a swimmer though such me- 
dia. While the equations of motions ([T]) remain invariant 
under Galilean transformation (up to a uniform pressure 
gradient), there is still an issue related to "permeability" 
of such matrix to a finite-size object propelled through 
it [i^l. The spacing between obstacles that compose the 
matrix places a constrain on the relevant dimension of the 
object moving through it without altering its structure. 
Obviously, for high volume fraction of the obstructions 
the effect of the matrix effect is most certainly local and 
not described by the effective mean resistance as in ([T|). 
The question to what extent the Brinkman's approxima- 
tion holds when the structure of the matrix is distorted in 
response to swimmer's movement, is yet to be answered. 
We can estimate the range of formal validity of the theory 
using some simple geometric arguments. For instance, 
for fibrous media {e.g. chain polymer aggregates in gels) 
where the fiber axes are completely randomly oriented 
the value of a is determined in a self-consistent manner 
as function of mean fiber content [37| . 

2 . - /^l 2 5 Ki{a,)\ 
V3 6 Ko(a*)J 

where a* = ab, (j) = \-k\? is the fiber volume fraction 
with A being the mean fiber length per unit volume and 
h stands for the fiber radius. Therefore, for sparse fiber 
matrix, </) = 0.001 (0.1% vol), we find that a* = 0.0304. 
Thus, for the typical value of ajh — 50 (e.5. the radius of 
the spirochete cell helix is ~ 0.05 /xm, the typical radius 
of polysaccharide chains in the extracellular domain is 
0.5 — 2 nm) the corresponding value of aa = ai,{a/b) = 
1.52 and the optimal velocity of propulsion is U/ujr ~ 
0.49 which is an improvement of « 40% comparative to 
the optimal speed of a helix propelled through viscous 
medium {U/ru> ~ 0.35). Assuming regular fiber spatial 
arrangements, one can estimate the average separation 
between fibers in the matrix A as 

A/6 « 2 r ^/^ - 1^ . 

For (j) = 0.001 this formula yields A/6 w 100, so that 
the mean spacing between fibers. A, is sufficiently large 
to assume that the propelled helix is not distorting the 
matrix microstructure. 



The situation is much more complicated when propul- 
sion in gel- forming semi-dilute solutions of some polymers 
(such as methylcellulose) are considered. The typical 
"mesh size" between polymer chain aggregates in these 
medium can be considerably higher than the above esti- 
mate based on uniform spatial arrangement, due to large 
concentration fiuctuations typical for semi-dilute solu- 
tions of fiexible and semi-fiexible polymers [3| • There- 
fore, the micro swimmer may navigate its way through 
polymer-lean region, taking advantage of the distributed 
hydrodynamic resistance due to clusters and/or aggre- 
gates in polymer-rich regions. The quantitative compar- 
ison of the present theory with the experimental results 
showing the significant increase of the propulsion veloc- 
ity in gel-forming polymer solutions is difficult due to the 
complex microstructure of the gels, the local interaction 
of the propelled object with the polymer network, elas- 
tic response of the distorted network, etc. It should be 
further emphasized that the present work is not a quan- 
titative description of propulsion in polymer gels, but it 
comes to qualitatively demonstrate the potential physical 
origin of propulsion enhancement due to screened hydro- 
dynamic interaction with the embedded network. 

The suggested numerical formalism of Sec. IIVBI 
can be readily extended to other cases of interest. If 
nearby obstacles are allowed to move in response to the 
propeller displacement when the certain stress threshold 
is exceeded, the rotating filament, for instance, could 
swim through arbitrary dense matrix of obstacles. Such 
approach can be applied for more realistic modeling 
of propulsion through viscoplastic materials exhibiting 
finite yield stress. A similar mechanism of locomotion 
enhancement is expected, as the swimmer is propelled 
through liquid-like region encapsulating the moving 
object, and it is "pushing" against the unyielded solid 
bilk material. In the limit of small volume fractions 
of obstructions the effective media approximation ([1]) 
should still hold, while in the other limiting case of freely 
suspended obstacles, the swimmer would experience 
hydrodynamic resistance controlled by the elevated 
effective viscosity, /icff, of the viscous suspension. Since 
zero-Reynolds-number propulsion through a viscous 
Newtonian liquid is purely geometric, propulsion ve- 
locity (for a prescribed swimming gait) is independent 
of the viscosity and no enhancement of locomotion is 
expected in this case. Adding elastic response of the 
distorted network of obstacles would allow more realistic 
modeling of propulsion in complex viscoelastic media 
such as viscous polymer gels. 
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